Anomalous Josephson current in superconducting topological insulator 
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We study the Josephson effect in superconducting topological insulator (STI) by referring to the 
pairing states of Cu-doped Bi2Se3. The current-phase relation of a junction between even-parity 
superconductors and STI shows the robust vanishment of the first order component, J(ip) ~ sin(2<£>), 
due to the mirror symmetry of the system. Furthermore, we clarify that the temperature dependence 
of critical current in junctions consisting of two STIs changes qualitatively depending on the relative 
spin-helicity of the two surface Andreev bound states. Detecting of these features qualifies as a 
distinct experimental evidence for the Majorana nature of the Andreev bound states. 

PACS numbers: 74.45.+C, 73.20.At, 03.65.Vf 



Introduction- Topological superconductors, which sup- 
port Majorana fermions as gapless Andreev bound states 
(ABS) on the surface or edge [lreli are a current hot topic 
in condensed matter physics. Recently, a new supercon- 
ductor, Cu-doped Bi2Sc3, dubbed as superconducting 
topological insulator (STI) has been discovered 0, [||. 
The presence of zero bias conductance peak in tunnel- 
ing spectroscopy supports the realization of topological 
superconducting state in this material 0, Eoj . The par- 
ent material Bi2Sc3 is a topological insulator, which is 
a band insulator with masslcss Dirac fermions on the 
surface 

@, EI El- Correspondingly, the STI supports 
surface Dirac fermions in the normal state. The Dirac 
fermions can be hybridized with ABSs in the supercon- 
ducting state, and as a result, the energy dispersion of 
Majorana fermions in the STI can be different from those 
in other topological superconductors The differ- 

ence in the energy dispersion enhances the surface den- 
sity of state (SDOS) at the Fermi energy, which affects 
on the tunneling conductance through the surface, re- 
markably [3 [H|. Moreover, several proximity effects 

and Josephson effect [20H30I ] of topological insu- 
lators have been reported. 

The purpose of this letter is to clarify Josephson effect 
of STI. Josephson effect is one of the most important 
quantum transport through surfaces in superconductiv- 
ity. Usually, they are proportional to sin^, where ip is 
the phase difference of pair potentials between two su- 
perconductors forming the junction. It is well known 
that 7r-phase shift [3l[ has been used for interferometers 
of Josephson junctions 32, |33|, w hich has established 
d-wave symmetry of cuprates 32h34|. In some d-wave 
superconductor junctions, the second-order component 
of the Josephson effect. sin2g?, can be dominant in the 
current-phase relation 35H38I]. Furthermore, Josephson 
junctions between spin-triplet superconductors and spin- 
singlet ones may show the second order behaviors as well 
39M42J]. The sin2y> behavior, which implies a double de- 



computers |43j. In the cases mentioned above, however, 
the second order behavior is easily lost due to the spin- 
orbit scattering at the interface. 

In this letter, we reveal that Josephson effect in STI 
are characterized by robust second order behaviors in 
the current phase relation. Based on the symmetry of 
the system, we first argue that the first-order component 
sin <p vanishes even in the presence of the spin-orbit scat- 
tering at the interface. Thus the leading term should 
be the second order. Then, we microscopically calculate 
the Josephson current for various junctions consisting of 
STI. For high T c cuprates, several interestin g fe atures of 
Josephson current have been found [3l|, [36M38L H ] due 
to ABSs f45M47j ] generated at the interface (surface). For 
STI, however, there has been no microscopic calculation 
of Josephson effect up to now. From the microscopic 
calculation, we establish the second order behavior of 
the current-phase relation, and clarify that the surface 
ABS enhances the Josephson current. Furthermore, we 
find an anomalous behavior of the Josephson current spe- 
cific to helical Majorana fermions in STI. In particular, 
it is shown that the magnitude of Josephson current be- 
tween two STIs is suppressed at low temperature when 
the spin-helicities of the helical Majorana fermions on 
the junction are mismatched. The anomalous tempera- 
ture dependence would provide a unique signal of surface 
helical Majorana fermions in STI. 

Symmetry-based argument- Let us start with a general 
argument based on symmetry of the system. Josephson 
current is generally decomposed into a series of different 
order, 



J((f) = ( J n sin ncp + I n cos nip) 



(1) 



n=l 



generate ground state of the junction, is of particular 
interest since it can be used for a quiet qubit of quantum 



Under the time-reversal, J(ip) goes to — J(— ip), thus in 
time-reversal symmetric junctions, J(ip) satisfies 

J(ip) = -J(-<p). (2) 

This implies that /„ = and the leading term is J(<p) ~ 
sin ip. 
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Now consider symmetry specific to STI. According to 
Ref.Q, there are four possible gap functions for STL 
Among them, however, only the two gap functions that 
belong to the A\ u and E u representations in the D 3( i 
point group are consistent with the tunneling conduc- 
tance experiment 0, [l5| . They both are invariant un- 
der time-reversal and they are odd under the inversion, 
which are responsible for topological superconductivity 
IE El S] • They are also odd under the mirror reflection 
with respect to the j/z-plane, (x, y, z) — > {—x, y, z). From 
the mirror symmetry, a structural transition in the en- 
ergy dispersion of Majorana fermions may occur, which 
enables us to explain the zero-bias peak of the tunneling 
conductance observed experimentally 
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We argue here that the mirror symmetry gives an ad- 
ditional constraint for J(ip) Consider a junction 
consisting of an s-wave superconductor and STI (s/STI). 
The interface of the junction is prepared so as the mirror 
plane of STI is perpendicular to it. Under the mirror 
reflection, the STI changes the sign of the gap function, 
while the s-wave superconductor does not. Consequently, 
one obtains the additional phase of 7r in the Josephson 
current, J{ip + Ti) 1 under the mirror reflection. Therefore 
the mirror symmetry implies 



J(ip) = J(^ + 7T). 



(3) 



This equation yields that Ji n +\ = ^2n+i = in Eq.([T]) 
and the first order term J(ip) ~ sirn^ vanishes. Thus the 
leading term is the second order J(ip) ~ sin 2<p in the 
s/STI junction. The 7r-periodicity of the Josephson cur- 
rent in Eq (0} is consistent with the effective Josephson 
coupling oc [(A*) 2 A| TI + h.c] between an s-wave gap A s 
and the STI one Agxi [8( . 

Our argument above implies that the second order be- 
havior in the s/STI junction is robust as far as the time- 
reversal and the mirror symmetries are preserved. In 
particular, it is not lost even when the spin-orbit scat- 
tering is present in the junction. This is different from 
previously known second order behaviors in Josephson 
current. For example, as is the case with STI, a chiral 
p-wave superconductor Si'2Ru04 is a topological super- 
conductor with chiral Majorana edge states 0,H3|> anc ^ 
the current phase relation between an s-wave and the 
chiral p-wave superconductor is proportional to sin2^, 
in the absence of spin-orbit interaction [HI], [52|. How- 
ever, the chiral p-wave superconductor is not odd under 
the mirror reflection, and also breaks the time-reversal 
invariance. Therefore, the second order behavior is frag- 
ile and actually the first-order term cos ip in the current- 
phase relation appears immediately in the presence of 
spin-orbit interaction [52j |. 

In the following, we also consider junctions consisting 
of d-wave superconductors and STI. As far as the inter- 
face is prepared so as the d-wave superconductor is even 
under the mirror reflection, we have a robust second order 
behavior of the Josephson current in a similar manner. 



Microscopic calculation- In the following we calculate 
Josephson current in STI microscopically. As regards 
the pairing symmetry of STI, we focus on the A\ u repre- 
sentation which is a full-gap topological superconductor, 
while qualitatively similar results are obtained for the 
nodal topological superconductor in the E u representa- 
tion. As mentioned above, the system has a diagonal 
mirror plane (x,y,z) —> (—x,y,z), and the gap function 
is odd under the mirror reflection. We consider various 
junctions, which consist of (a) s-wave superconductor, 
normal metal (N), and STI (s/STI), (b) d yz -wave super- 
conductor, N, and STI (d yz /STT), (c) d y 2 _ z 2-w&\e super- 
conductor, N, and STI (d y 2_ z 2 /STI), and (d) STL N and 
STI (STI/STI), respectively. The orientation of the su- 
perconductors in the junctions are chosen as illustrated 
in Figs. [T] and El 

The calculation is based on the three-dimensional 
square lattice model. We put s, d yzi d y 2_ z 2- 
wave superconductor, or STI in the left re- 
gion (I < z < Al), N in the center region 
(A L + 1 < z < A L + A c ), and STI in the right 
region (A L + N c + 1 < z < A L + A c + A R ), respectively. 
The on-site energy c^tACn and the nearest neighbor hop- 
ping along the z-direction c' n tACn+i with A = s, d, STI 
are given by e s = [2t(cosk x + cosk y ) — f^ s ]T Zl e.d — 
[2iy cos fc x + 2£j_cosfc y - Mrf] T z,£STi = {[^o + 2mi + 
4to2 — 2to2 (cos k x + cos ky )] CF x +va z {sy sin k x — s x sin k y ) — 
MSTij-Tz, t s = tr z ,t d = t||T z ,tsTi = (mi<j x — i^o-y) r 2 , 
where a, s, and r are the Pauli matrices in the orbital, 
spin, and Nambu spaces, respectively. In s/STI (d/STI) 
junction, the on-site energy cn and the hopping t-^ in N 
are chosen to be the same values as those in the s(d)- 



wave superconductor: 



s(d)j *n 



t, 



(<*)• 



In STI/STI 



junction, we set 6n = e s in N. For the cases (a) to (c), 
the hopping between the s(d)-wave superconductor and 
N is assumed to be the same as that in the s(d)-wave 
superconductor. On the other hand, in STI/STI junc- 
tion, the transfer at the left interface between STI and 
N is given by H Th = J^as ci as thaT z c n+ls + h.c, 
with n = Al, and at the right interface, the 
transfer is H TR = ^] (Ts c\ ls t Ka T z c n+las + h.c, with 
n = Al + Nr. Pair potentials of s (A s ), d yz (A yz ), 
d y 2_ z 2 (A, y 2_ z 2) -wave superconductors, and STI 



(Asti) are given by A s 



E 



-1 AsTjcCn, A yz 

h.C, Ay2_j.2 



S n il 1 c n~i&d Sin k y T x C n+1 
Er^=l cK- A d COS k y )T x C n + (j^n^l 1 c X^fc n+1 + h.C. J , 

and 

Asti = c^Asno-yS^Tx cos tp - T y sinip)c n . 

n=N Ij +N c + l 

(4) 

The temperature dependence of the pair potentials is 
given as that in the weak coupling limit: A(T) = 
Atanh(1.74^/T c /T - I), with A being the pair poten- 
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(a) s /STI 




FIG 

d, 



1. Geometries of (a): s/STI, (b): d^/STI, and (c): 
_ Z 2/STI junctions. We assume thin normal layers (black 



shaded region) at the center of the junctions. The Josephson 
current flows along the z-direction. 
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FIG. 2. The current-phase relations of STI junctions, (a) 
s/STI junction with non-topological STI. Here a conventional 
s-wave pairing is realized in STI. (b) s/STI junction, (c) 
d V z /STI junction, (d) d y i_ z i /STI junction. In (b)-(c), topo- 
logical STI is realized. The pairing in the topological STI is 
given by Eq. Q. We use the parameters taken from Ref. [l3| , 
as follows. a a = 0, tru = tn2 = —0.1, t = —1, A ST i = 0.1, 
A s = 0.2, m = -0.7, mi = 0.5, m 2 = 1.5, ^sti = 0.9, 
v z = 1, v = 1.5, N L = N R = 20, N c = 2, and T = 0. The 
parameters of the two d-wave superconductors are taken as 
follows. U = -1, t± = -0.1, fi d = -0.5, A d = 0.5. 



tial at the zero temperature: A = 1.76T C [28|, [53|. Un- 
der these settings, the Josephson current is calculated 
h y J = 2^ 3 'Efc||(<4(*!||)*NC„-(-i(fc||)> + c.c, where iV| 
is the number of the unit cells in the xy-planc, (• • • ) 
is the thermal average, and n is a site in N, that is 
n e [N h + l,N L + N c -l]. 

The obtained current-phase relations in s/STI, 
dy Z /STl, and d y 2_ z 2 /STI junctions are summarized in 
Fig. [2] As a reference, we also plot the current phase 
relation of non-topological s/STI junctions where con- 
ventional s-wave pairings are realized in STI. It is found 
that the Josephson current shows the second order behav- 
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FIG. 3. Josephson current in s/STI junction under in-plane 
Zeeman field, (a) Geometry of the system, (b) ^-dependence 
of Fourier components of the Josephson current. The sin ip 
and cos 2yp-components are identically zero. See text, (c) 6h- 
dependence of the local minimum values (at cp — ±7r/2) of 
the free energy f(ip). The strength of Zeeman field is chosen 
as V a sti = 0.1. 



ior J ((f) <~ sin2^j, except for the non-topological s/STI 
one. J (if) ~ sin<y9 is satisfied for the non-topological 
s/STI one. In the non-topological case, the gap function 
in STI is even under the mirror reflection, thus cq.(|3|) 
does not hold. Consequently, J(<p) has a conventional 
sinusoidal dependence. On the other hand, J ((f) in the 
topological STI cases has a sin 2^ dependence. It is re- 
markable that the magnitude of the Josephson current 
in dy Z /STl junction is much larger than that in s/STI or 
dj,2_.j2/STI junction. This is because there exists a flat 
band ABS of d yz -wa,ve superconductivity |45j on the in- 
terface of dy Z /STl junction, which enhances the Joseph- 
son current at low temperature [3(1, |4(| 5^ . The sin 2tp 
behavior (or the vanishment of the first orders) is a dis- 
tinct character of the topological superconductivity in 
STI. 

By breaking the mirror and time reversal symmetries, 
the first order terms in the current phase relation can be 
induced. Indeed, under an in-planc Zeeman field Hz = 
Y^n c n h ' sc " witn h x = 0,h y = hsin9h,h z = hcos9 h , 
Eqs.((2]) and ((3]) change as J(<p,0h) = — J(— f, O^ + tt) and 
J(ip, Oh) = J(<p + 7T, Oh + tt), respectively, with J(ip, Oh) 
the Josephson current under the Zeeman field. From 
them, we have J{f, Oh) = — J(— </? + 7r, Oh) as a constraint. 
Therefore, as confirmed numerically for s/STI junction 
under a Zeeman field illustrated in Fig. 03a), the first 
order component cos ip can appear in the current-phase 
relation, while siny and cos2ip can not [Fig. EJb)]. 

Now discuss the free energy of Josephson junctions of 
STI. The free energy / per unit area of the Josephson 
junction is given by f(ip) = h/(2e) J Q dtp' J(<p')- There- 
fore, the second order behavior J((p) ~ sin2<p implies 
doubly degenerated ground states in the STI junctions. 
The minima of f(<p) arc located at tp = ±7r/2. As dis- 
cussed in the above, by applying the Zeeman field, the de- 
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generacy can be lifted. Figure [3jc) shows 6^-dependcncc 
of f(f) at (p = ±7r/2, where the free energy /(</?) takes 
local minima. This figure illustrates that one may choose 
either of the ground states by applying a suitable Zeeman 
field. 

Spin-helicity dependent Josephson effect- Finally, we 
clarify Josephson effect intrinsic to helical Majorana 
fermions in STI. As mentioned above, the STI supports 
helical Majorana fermions on its surface. The effective 
Hamiltonian of the surface helical Majorana fermion on 
the interface of junctions (the xy plane) is represented as 



(a) STI(+ 



(b) STI(- 



H SU if (fe|| ) — w surf {^xSy kyS x ), 



(5) 



near fc|| = 0. From the Hamiltonian, one can obtain 
the linear dispersion of Majorana cone, E SUT f = ±v snl f k\\ 
with fcii = |fc|||. The spin and the momentum are locked 
on the cone so as the Majorana fermions have a definite 
eigenvalue h s of the spin-helicity (fcii x s)/ku in low en- 
ergy. From eq.([5]), the spin-helicity is h s = sgn(w sur f) 
[h s = — sgn(v sur f )] for the upper (lower) cone. A unique 
character of Majorana fermions in STI is that the spin- 
hclicity can change depending on the chemical potential 
of the system. As shown in Fig. SI when the chemical po- 
tential /xsti increases, the spin-helicity of the upper cone 
near fc|| = changes from h s = — to h s = + at a critical 
value /^| TI . We also find that when /xsti < Mstd the en- 
ergy dispersion of the Majorana fermion is not a simple 
cone, but a rather complicated caldera-shaped one. As a 
result, in addition to the cone with h s = — near fey =0, 
there is another branch with h s = + in the spectrum, as 
illustrated in FigJH Referring to the spin-helicity of the 
upper cone near fc|| = 0, we denote STI with a simple 
cone and that with a complicated caldera shaped one as 
STI(+) and STI(-), respectively. 

Now let us study the spin-hclicity dependence of 
Josephson current. For this purpose, we consider 
STI/STI junction illustrated in Fig. EJ By tun- 
ing the chemical potentials in the left and right STIs, 
we can change the spin-helicity of the Majorana cone 
in each STI. We calculate the Josephson current of 
STI/STI junctions with three different combinations 
of spin-helicity, STI(+)/STI(+), STI(+)/STI(-) and 
STI(— )/STI(— ), respectively. In these junctions, the 
current-phase relation of J(<p) is rather conventional, 
however, the temperature dependence becomes anoma- 
lous. Figure [5] shows the temperature dependence of 
the maximum Josephson current max v ( J) for these junc- 
tions. We find that, as compared to the junctions 
with the same spin-helicity, i.e. STI(+)/STI(+) and 
STI(— )/STI(— ) junctions, the Josephson current with 
the mismatched spin-helicity [STI(+)/STI(— ) junction] 
is strongly suppressed. In particular, we find that the 
Josephson current of the latter junction decreases at low 
temperature while those of the former increases. Here 
we note that the suppression in the mismatched case oc- 
curs below T/Tcsti, ~ 0.07, which exactly corresponds 




0.1 0.2 0.3 0.4 0.5 0.1 0.2 0.3 0.4 0.5 0.6 

k x k x 

FIG. 4. Energy spectra of the surface states in STIs. (a) The 
positive energy part of the surface state for usti = v z in the 
region of k y = and k x > 0. Inset shows the corresponding 
overall spectrum of the h s — + branch of the surface state, 
(b) The positive energy part of the surface state for ^sti = 
0.9v z in the region of k y = and k x > 0. Inset shows the 
corresponding overall spectrum of the h 3 = — branch of the 
surface state, which is discussed in Ref. flal. 
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FIG. 5. Temperature dependences of Josephson current J(T) 
in STI/STI junctions. Geometry of the system is illustrated 
in the inset. Here, STI(+)/STI(+) [STI(-)/STI(-)] junc- 
tion is prepared by choosing the chemical potential as /ig TI = 
Msti = v z (msti = Msti = 0.9^), where Usti and ^g TI are 
the chemical potentials in the left and right STIs respectively. 
The spin-helicity mismatched case of STI(+)/STI(— ) junc- 
tion is constructed by putting /xg TI = v z and /isri = 0.9w z . 
The values of hoppings between N and STIs are chosen as 
thi = £l2 = £ri = £r,2 = O.lt. 



to the energy \E/A\ < 0.07 where the h s = — branch 
appears in STI(-) [See Fig. H[b)]. This implies that 
the suppression of the Josephson current occurs due to 
the mismatch of the spin-hclicity in the junction. Wc 
also find that the Josephson current in STI(— )/STI(— ) 
is much more enhanced at low temperature than that 
in STI(+)/STI(+) because a twisted energy spectrum of 
the caldera cone has many low-lying states which con- 
tribute to Josephson current. We would like to mention 
here that the spin-helicity dependence of the Josephson 
current is different from that of two-dimensional helical 
superconductors since the Josephson current in two di- 
mensional helical superconductors is always enhanced at 
low temperature, independent of the spin-hclicity [54 1. 
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Discussions- The anomalous Joscphson effect reported 
in this letter are accessible experimentally: First, the 
vanishment of sin ip- and cos pterins in the current-phase 
relation in s/STI and d/STI junctions is detectable by the 
Shapiro steps at the bias voltage given by V = nhw/{_Ae), 
where n is an integer and ui is the frequency of the mi- 
crowave. The DC-SQUIDs with these junctions are also 
very sensitive to the second order behavior sin 2(p of the 
current-phase relation. Moreover, the anomalous tem- 
perature dependence of STI/STI junctions is easily mea- 
surable in experiments. The anomalous temperature de- 
pendence is of particular interest since it is a direct ex- 
perimental signal of the spin- lock nature of surface helical 
Majorana fcrmions, and thus a direct experimental evi- 
dence of topological superconductivity of STL 

To summarize, we have shown that Josephson current 
between even-parity superconductors and STI can be ro- 
bustly proportional to sin2ip as far as the mirror sym- 
metry of the system is preserved. Therefore, Cu^T^Scs 
is a good platform to realize a Josephson junction with 
sin 2<ys-current-phase relation and two-fold degenerated 
ground states, which can be utilized for quiet qubit in 
quantum computing. We have also clarified that the 
magnitude of the Josephson current in STI/STI junction 
is influenced critically by the relative spin-helicity be- 
tween the surface helical Majorana fermions of the two 
STIs. Although the existence of s urface gapless states 
has been experimentally confirmed in Cu^I^Sea by tun- 
neling spectroscopy 0], which is consistent with topo- 
logical superconductivity in Cu K Bi 2 Se3 [3 IH, HBJ] j their 
Majorana nature has not yet. The spin-helicity depen- 
dent Joscphson current could be an experimental signal 
that is directly related to the Majorana fermions on the 
surface of STI. 

This work is supported by the "Topological Quantum 
Phenomena" (No. 22103005) Grant-in Aid for Scientific 
Research on Innovative Areas from the Ministry of Edu- 
cation, Culture, Sports, Science and Technology (MEXT) 
of Japan. 
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